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ABSTRACT 

We present low-resolution MMT Hectospec spectroscopy of 594 candidate Monoceros stream mem- 
ber stars. Based on strong color-magnitude diagram overdensities, we targeted three fields within the 
stream's footprint, with 178° < I < 203° and —25° < b < 25°. By comparing the measured iron 
abundances with those expected from smooth Galactic components alone, we measure, for the first 
time, the spectroscopic metallicity distribution function for Monoceros. We find the stream to be 
chemically distinct from both the thick disk and halo, with [Fe/H] = — 1, and do not detect a trend in 
the stream's metallicity with Galactic longitude. Passing from 6=+25° to &=— 25°, the median Mono- 
ceros metallicity trends upward by 0.1 dex, though uncertainties in modeling sample contamination 
by the disk and halo make this a marginal detection. In each field, we find Monoceros to have an 
intrinsic [Fe/H] dispersion of 0.10-0.22 dex. From the CallK line, we measure [Ca/Fe] for a subsam- 
ple of metal poor program stars with —1.1 < [Fe/H] < —0.5. In two of three fields, we find calcium 
deficiencies qualitatively similar to previously reported [Ti/Fe] underabundances in Monoceros and 
the Sagittarius tidal stream. Further, using 90 spectra of thick disk stars in the Monoceros pointings 
with b ~ ±25°, we detect a 0.22 dex north/south metallicity asymmetry coincident with known stellar 
density asymmetry at Rgc ~ 12 kpc and \Z\ « 1.7 kpc. Our median Monoceros [Fe/H] = —1.0 and 
its relatively low dispersion naturally fit the expectation for an appropriately luminous My ~ —13 
dwarf galaxy progenitor. 

Subject headings: Galaxy: evolution, Galaxy: Stellar Content, Galaxy: structure, galaxies: dwarf, 
galaxies: interactions, stars: abundances 



1. INTRODUCTION 

The Monocero s stream (Ncwbcr g et al.l 120021 : 
lYannv et all [20031 hereafter Y03) comprises a ring- 
like stellar overdensity in the plane of the Galactic 
disk. Between 110° < I < 250°, the kinematically cold 
(Y03) stream covers —30° < b < 30° at galactocentric 
radius R~17-19 kpc with radial thickness AR ~ 4 kpc 
(| Juric et al.ll2008l ). 

The true nature of the structure is still uncertain, with 
the leadin g scenarios being: a) remnant of an accre- 
tion event (IPenarrubia et aLl|2005|) . with the Canis Major 
dwarf galaxy |Martin et al.l 120041 ) as a possible progen- 
itor, b) disturbance due to a high-eccentricity fly by en- 
coun ter (lYounger et al.l I2Q08T) . and c) disk flare or warp 
(e.g- JMomanv et al.ll2006D 

Discerning between these scenarios, especially the 
first two, may have important theoretical consequences 
as mergers with orbit s and mass ratios implied by 
IPenarrubia et all (I2005P are deemed unli kely in current 
ACDM models (e.g.. lYounger et al.ll200l . 

In characterizing any stellar population, metallicity 
is considered a fundamental parameter. Unfortunately, 
the literature regarding the chemical composition of the 
Monoceros stream is sparse and seemingly inconsistent. 
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For example, Y03 estimate [Fe/H] 1.6 at I = 198°, 

b = 27° using the Ca II K line a n d colo rs of turnoff stars. 
On the other end, iCrane et al.l (|2003[ ) (henceforth C03) 
measure a mean of —0.4 ± 0.3 using M giants, sampled 
mainly from observations of the northern component of 
the s tream. Mos t recen tly, using photometric metallic- 
ities, llvezic et al.l (|2008l) (hereafter 108) find the MSTO 
stars in the northern anticenter portion of the stream to 
have a median metallicity of [Fe/H] = —0.95, with intrin- 
sic RMS scatter of only 0.15 dex. Hypotheses invoked to 
reconcile t hese values range from m ultiple stellar popula- 
tions (e.g.. IPenarrubia et al.l 120051 ) . to systematic errors 
in calibrations (C03), and to questioning whether the 
northe rn and southern part of the ring share the same 
origin (jConn et al.l 120071 ). Also, except for the 108 pho- 
tometric study, the full metallicity distribution function 
(MDF) of stars in the Monoceros stream has actually 
never been measured from spectroscopy. 

To better characterize the chemical composition of 
Monoceros, we have obtained spectra of ~600 candidate 
member stars. This large sample size allows us to probe 
the Monoceros MDF with good statistics. Also, by point- 
ing in multiple directions along the stream, we can search 
for metallicity gradients. Specifically, we set out to char- 
acterize the [Fe/H] MDF of Monoceros in each of three 
fields with 177.9° < I < 202.9° and -24.2° < b < 24.4° 
(see Table Q] for field details and naming conventions) . 
Additionally, we aim to gauge the a-enhancement of 
Monoceros from a [Ca/Fe] measurement. This allows 
us to address the origin of Monoceros by comparing its 
[Ca/Fe] to the differing a-abundance patterns found in 
the Milky Way versus dwarf spheroidal (dSph) galaxies. 

2. MMT OBSERVATIONS 
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TABLE 1 

Summary of MMT Observations/Targets 



Sample 


I 


6 


a 


S 


Exposures (s) 


N stars 


(S/N) med 


^med 


Z|med a (kpc) 














NORTH 












Monoceros 


187.245° 


24.401° 


115.741° 


32.662° 


2x7200 


199 


49 


19.55 


3.871 


0.12 


Thick Disk 


187.245° 


24.401° 


115.741° 


32.662° 


2x7200 


30 


123 


17.62 


1.603 


0.12 












SOUTH 












Monoceros 


202.884° 


-24.183° 


76.908° 


-2.587° 


2x7200 


146 


17 


19.83 


4.911 


0.30 


Thick Disk 


202.884° 


-24.183° 


76.908° 


-2.587° 


2x7200 


60 


137 


17.44 


1.788 


0.30 












NORTH18 












Monoceros 


177.879° 


18.055° 


104.834° 


38.852° 


2x7200,5400 


249 


84 


19.08 


2.688 


0.29 












M13 












Calibration 


58.999° 


40.917° 


250.416° 


36.454° 


3600 


189 


39 


18.34 


5.028 


0.05 



a using the photometric parallax relation of 108 equations (A6) and 

2.1. Sample Selection 

Recognized as a stellar overdensity between b = ±30° 
for 110° < I < 250°, Monoceros appears most strikingly 
in the number counts of MSTO F/G dwarfs at galacto- 
centric distances R^15-20 kpc. Near the Galactic an- 
ticenter, these stars correspond to apparent magnitude 
r r; 18.5-20.5 and color g — r w 0.25-0.4, with ugriz re- 
ferring to SDSS photometric passbands. By examining 
(g, g — r) color-magnitude diagrams (CMDs) of objects 
classified by the SDSS (DR7) as stars within the joint 
Monoccros/SDSS footprint, we identified many degree- 
scale pointings in which the Monoceros MSTO was ob- 
viously visible above the Galactic background (for back- 
ground estimation details see 33. 6[) . In the best cases, 
the relevant CMD region of such pointings exhibited a 
factor of ~2 overdensity, with total projected density 
^200—300 stars/deg 2 . This field size and target density 
matches that of MMT/Hectospec, which, with its large 
6.5m mirror size, can obtain at r = 20 S/N sufficient to 
measure [Fe/H] with integration times of ~4 hours. 

Hoping to gauge both I and north/south gradients in 
the stream's metallicity, six such fields were initially tar- 
geted for MMT/Hectospec spectroscopy, covering 150° < 
I < 230° in the Galactic north and sampling four Galac- 
tic latitudes with —30° < b < 30°. In several of these 
fields, a relatively small subsample of much brighter thick 
disk targets were chosen from a similar range of MSTO 
g — r colors and with 17 < r < 18. 

2.2. Observations and Data Reduction 

The observations were carried out in MMT queue- 
scheduled multi-object observing mode. Due to a combi- 
nation of factors including weather, only three fields were 
observed, spanning a limited range of (I, b), though still 
including both Galactic north and south. Table Q] sum- 
marizes the targets/observations, and Figure [1] shows the 
CMD of each field, indicating targeted regions. 

Hectospec can observe up to 300 objects at once. In 
general, we typically devoted ^50 fibers to measuring the 
sky, with the remainder targeting Monoceros candidates 
subject to the constraints of the instrument (e.g., the 
availability of the targets in the desired magnitude range, 
and limits on how closely spaced the fibers could get). 
Any remaining fibers were set to target the thick disk 
stars, to investigate the thick disk asymmetry present 
towards the Galactic anti-center (see 35. 2ft . We used the 
270 gpm grating yielding a resolving power of R ~ 1,000 
near 5000A. The low resolution warrants full wavelength 



(A7) 

coverage of 3700-9150A. 

The data were reduced with customized IRAF proce- 
dures as part of t he SAO Telescope Data Center service 
(jMink et al.ll2007n . The S/N ratio of the final spectra is 
typically ^50-80 per pixel at ~ 5200A for the Monoceros 
stars and ^135 for the brighter thick disk stars. We used 
the IRAF/rvsao routine xcsao to measure the radial ve- 
locity (v rac i) of each star from cross-correlation with tem- 
plate stellar spectra on a per-exposure basis, correcting 
for heliocentric velocity before stacking the exposures of 
each object. Typical errors on v rac i ranged from 10-20 
km s _1 , varying based on S/N. 

2.3. Stellar Parameters 

Abundance extraction requires estimates of several 
stellar parameters of each object, namely: the surface 
gravity logg, temperature T e ff, microturbulcnce vt, and 
metallicity. As our targets were selected from SDSS 
imaging, each star has ugriz photometry available. Us- 
ing 108 equation (3), which gives Tcff as a function of 
g — r, we determined effective temperatures for all pro- 
gram stars. The statistical error per star due to photo- 
metric errors is ^80K for thick disk targets and ^110K 
per star for the fainter Monoceros candidates. T c g gen- 
erally ranged from 5800-6300K within each sample of 
Monoceros/thick disk stars. 

Our targets were selected to be MSTO F/G dwarfs. 
Based on the stellar parameters of 46 dwarfs from 
iFulbrighti pOOOh with -2 < [Fc/H] < 0.0 and 5800 
K< T eff <6300 K, we adopt v t = 1.25 km s" 1 and 
log<? = 4.20 for all Monoceros candidates and thick disk 
stars. In 33.5l we discuss the impact of these assumptions. 

Of course, we wish to derive rather than assume metal- 
licities. To extract abundances, we first determine [Fe/H] 
from iron absorption features assuming [X/Fe]=0 for all 
X, choosing features where this assumption has negligi- 
ble effect (see Figure |2(a)| and g33J). The [Fe/H] value 



of each object is then used as the metallicity input for 
subsequent Ca abundance measurements. 

3. IRON ABUNDANCE OF THE MONOCEROS STREAM 

To infer iron abundances from our spectra, we fo- 
cus on the iron line complex near A « 5270A. The 
absorption equivalent width (EW) in the region with 
5267A < A < 5287A is dominated by Fe, mainly Fe I, but 
with some Fe II. The primary Fe lines of interest and ma- 
jor conta minat ing lines due to other elements are labeled 
in Figure 2(a) All abundances, including those drawn 
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Fig. 1.— SDSS DR7 CMDs of target fields NORTH (a), SOUTH (b) and NORTH18 (c). All three panels have the same dynamic range 
(0-30 stars/pixel, with the pixel sizes being 0.025 mag and 0.1 mag in g — r and r, respectively). The black dots represe nt Monoceros 
candi dates from which the spectroscopic sample we analyzed was drawn. The black line shows an 8 Gyr, Z = 0.0039 isochrone (Mari go et al.l 
2008) intended only to guide the eye along the Monoceros MSTO. 



from the literature for calibration, wer e standardized 
with r espect to the current solar values of lAsplund et al.l 

(pool . 



An example best-fit result for typical S/N and iron ab- 
sorption strength is plotted in Figure 2(a) The line list 
employed in our syntheses is provided in Table El 



3.1. Continuum Normalization 

Given our relatively low resolution and moderate S/N, 
continuum normalization is a sensitive matter. We base 
our continuum fitting procedure on that us ed by the 
SEGU E Stellar Parameter Pipeline (SSPP, iLee et al.1 
2008) to measure [Fe/H], since our spectra have simi- 
lar S/N and resolution to those for which their pipeline 
was designed. In particular, we iteratively fit the region 
of the spectrum blueward of 5800A with a ninth order 
polynomial, masking the strong Balmer lines. The SSPP 
additionally computes abundances with local continuum 
estimates based on sidebands. However, our resolution is 
a factor of ^2 lower than that of SSPP spectra, making 
it unfeasible to isolate sideband regions which attain the 
continuum and are sufficiently wide to be robust. 

3.2. [Fe/H] Fitting Technique 

We take the calculated T e ff and assumed logg, v t for 
each star and compare the observed spectrum to a library 
of synthetic spectra generated based on these fixed stellar 
parameters, but with varying [Fe/H]. In particular, syn- 
thetic, continuum normalized spectr a were generat ed us- 
ing MOOG (the latest version, 2010: lSnedenlll973[ ). with 
one-dimensional no-overshoot model atmosphe res (as- 
sumin g local thermodynamic equilibrium) from iKuruczl 
( jl993f ). This version of MOOG accounts for the fact 
that Rayleigh scattering becomes an important source 
of continuum opacity at short wavelengths blueward of 
4500 A (|Sobeck et al.ll2011l) . For each [Fe/H] value with 
—3.5 < [Fe/H] < 0.5 in gradations of 0.01 dex, we com- 
pute the x 2 between the synthetic spectrum and the ob- 
served, normalized spectrum over the Fe absorption re- 
gion. We assign each star the [Fe/H] corresponding to the 
minimum % 2 value thus obtained. The \ 2 is computed 
with the variance taken from Poisson statistics to be 
a 2 = (Nf g+R 2 ), where N is the counts per pixel in ADU, 
R = 2.8 ADU is the readnoise and g = 1 is the gain. We 
compute the \ 2 within the range 5267A < A < 5287A, 
a region determined to optimize the Fe absorption rela- 
tive to absorption by other elements and random noise. 



3.3. Iron Abundance Calibration 

To test our [Fe/H] fitting procedure, we apply the 
above methodology to 189 MMT/Hcctospec spectra of 
globular cluster stars in M13 with comparable S/N and 
SDSS photometry. We again estimate T e g from g — r 
and adopt v t = 1.5 km s _1 , logg = 4.25 based on 
iBrilev fc Cohenl (120011) Table 2. Ml 3 is kn own to have 
[Fe/H] = -1.55 (jCohen fc Melendezl 120051 ) . and we re- 
cover a median [Fe/H] of —1.51. We thus apply a —0.04 
dex correction factor to all of our [Fe/H] measurements. 

3.4. Rejecting Poorly Determined [Fe/H] 

Based on visual inspection of the faintest and appar- 
ently most metal poor objects in our Monoceros candi- 
date sample, it became clear that several of these lacked 
sufficient iron EW to measure [Fe/H] in the presence 
of statistical noise and minor continuum placement er- 
rors. To remove these outliers, we computed ctewi the 
la error in total Fe fitting region EW due to Pois- 
son noise and continuum placement errors. The con- 
tinuum error was estimated by fitting a polynomial two 
orders lower than the default continuum, and generally 
amounted to <1% continuum level difference. We then 
rejected all stars for which the S/N and EW were suf- 
ficiently low that (Jew translated to >0.5 dex [Fe/H] 
uncertainty. For Monoceros targets in fields NORTH, 
SOUTH, NORTH18 the percentage of stars rejected was 
7.5%, 6.8% and 0.4% respectively. All rejected stars had 
estimated [Fe/H] < -2.05. No thick disk [Fe/H] esti- 
mates were rejected, as these stars generally have higher 
S/N and metallicity than the Monoceros candidates. 

In performing maximum likelihood fits to our [Fe/H] 
distributions (see £ 13 . 6[) we also identified, among Mono- 
ceros candidates, three [Fe/H] outliers at [Fe/H] > 0. 
These objects also displayed unreasonably low [Ca/Fe] < 
— 1, again suggesting a dramatic overestimate of [Fe/H]; 
these three objects were also excluded from the sample. 

3.5. Iron Abundance Uncertainties 

There are several uncertainties on our derived metal- 
licities. First there is statistical scatter from Poisson 
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Fig. 2. — (a) An example of the [Fe/H] fitting region for typical [Fe/H]. The blue line with errorbars is the observed data, the green line 
is the best-fit synthetic spectrum, the red dashed line is the continuum fit, and the purple dashed lines show the fitting region boundary. 
The upper/lower green dashed lines show synthetic spectra with [Fe/H] lowered/raised by 0.25 dex relative to the best fit. The strong Fel 
line at 5269. 5A, plus the next 9 most important lines by EW are marked with black vertical lines, (b) Same for CallK at 3933. 7k, plus 
next ten most important lines by EW, all of which are Fef features. The upper/lower green dashed lines correspond to offsets in [Ca/Fe] 
of =f0.25 dex relative to the best fitting [Ca/Fe], keeping [Fe/H] fixed. 



noise, which, for typical S/N is a function of [Fe/H]. We 
gauge this error with Monte Carlo generation of spectra 
to which we have added random Poisson noise. For typ- 
ical [Fe/H] = -1, we find an RMS scatter of 0.11 dex. 
At typical halo [Fe/H] = —1.6 we find a scatter of 0.25 
dex, while for typical thick disk [Fe/H] = —0.7 we find 
an RMS scatter of only 0.05 dex. With sample sizes of 
100—250 stars per field, this statistical error is thus not 
a primary concern. 

Another source of uncertainty stems from our assump- 
tions about stellar parameters. We find the effects of 
reasonably large systematic over/underestimates in each 
of \ogg (±0.3 dex), T cS (±100 K), and v t (±0.3 dex), to 
be T0.05 dex, ±0.07 dex, ±0.07 dex in [Fe/H] respec- 
tively. Summing these effects in quadrature yields an 
uncertainty of ±0.11 dex on our derived metallicities. 

A further possible source of bias in [Fe/H] is contam- 
ination of the iron abundance fitting region by lin es of 
other elements, most notably Ca (see Figure 2(a)). In 
determining [Fe/H], we generated synthetic spectra as- 
suming that [X/Fe]=0 for all X. While [X/Fc] is not 
known for Monoceros, the most relevant parameter of 
[Ca/Fe] for the Galactic disk/halo tends to lie between 
0.0 and 0.4. We address the resulting [Fe/H] bias by us- 
ing MOOG to simulate populations of stars with [Ca/Fe] 
following the known Galactic trend with [Fe/H] (see 
§4.ip . adding appropriate noise to the synthetic spec- 
tra before feeding them back into our [Fc/H] fitter. The 
resulting bias is very small, <0.02 dex. 

A final important source of uncertainty is the con- 
tinuum normalization procedure. With our low reso- 
lution, the blending of lines means the continuum is 
rarely achieved, and could potentially be systematically 
over/underestimated throughout the sample, thereby 
raising/lowering all metallicities. To test the system- 
atic effect due to differing continuum normalizations, we 
again re- fit abundances with new continua parameterized 
by a polynomial two orders lower than the default con- 
tinuum. We find in general a resulting systematic offset 
of 0.07 dex in [Fe/H]. 



We assign a final systematic uncertainty of 0.13 dex to 
our [Fe/H] distributions, adding the errors from contin- 
uum fitting and stellar parameters in quadrature. This 
value is somewhat conservative, as much of this uncer- 
tainty (particularly the continuum normalization) has 
been calibrated out with the M13 data. We find, through 
simulations similar to those previously described, that 
the observed [Fe/H] scatter is well described by a com- 
bination of both Poisson noise and statistical scatter of 
the underlying stellar parameters, again at the level of 
±100K in T c ff, ±0.3 dex in both v t and \ogg. Thus, 
summing these sources of scatter in quadrature, we find 
for typical S/N and [Fe/H] = — 1 a single-star statistical 
error of 0.16 dex. 

3.6. Decomposing the [Fe/H] Distributions 

To characterize the Monoceros MDF, we must first 
decompose the distribution of reliable [Fe/H] estimates 
in each field into Galactic components (halo, disk) 
and Monoceros members. To estimate the expected 
background in both halo and disk stars, we simulate 
CMDs for each po inting with both the Besancon model 
(jRobin et al.ll2003f) and qalfa s t, usin g the best fit Galac- 
tic parameters of ljuric et al.l (|2008f ). All simulations in- 
cluded realistic photometry errors and were extinction 
free, for com parison to observed SP SS CMDs dereddened 
according to iSchlegel et al.l (|1998ft . henceforth SFD. We 
created a matched filter in target density for the stars 
with reliable [Fe/H] in each field, then integrated against 
the observed, Besangon and galfast star counts to deter- 
mine the expected number of contaminating stars from 
the disk and halo. For each field, overall CMD normal- 
ization constants close to unity were applied by requir- 
ing agreement between the total number of stars redder 
than g — r = 0.8 and falling within the targeted range 
of r values. This analysis suggests ^65% of Monoceros 
candidates with well-measured [Fe/H] are true stream 
members. In no case do we expect any thin disk stars, 
as the Monoceros targets have |Z| > 2.5 kpc. 

Figure [3] shows a decomposition of the observed dis- 
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tributions of reliable [Fe/H] for each field. In princi- 
ple, there are nine free parameters: the normalization 
of Monoceros, thick disk, and halo components, their 
median metallicities and their dispersions. The follow- 
ing subsections describe our assumptions for the thick 
disk and halo distributions, and procedure for fitting the 
Monoceros component. 

3.6.1. Thick Disk Parameters 

In each field, the thick disk normalization was fixed 
by the CMD analysis to be the average of the predic- 
tions from the galfast and Besangon simulations. In fields 
NORTH, SOUTH we fix the median metallicity to that 
of the thick disk tar gets in these fields, plus a small —0.02 
dex correction from lBond et al.l (|2010D equation (A2) to 
account for the vertical disk metallicity gradient (our 
thick disk targets are at lower \Z\ than the more dis- 
tant thick disk stars which contaminate the Monoceros 
sample). For field NORTH18, in which we did not tar- 
get any thick disk stars, the median [Fe/H] used for field 
NORTH was adopted. A 0.35 dex intrinsic dispersion cal- 
cu lated from th e kinem atically selected thick disk stars 
m iBensbv et all (f2007h was adopted for the thick disk in 
all fields. 

3.6.2. Halo Parameters 

In fields NORTH, NORTH18, the halo normalization 
was fixed to the average of the galfast and Besangon 
predictions. A cut on \v ra d\ > 150 km s _1 , yielding a 
subsample which should be dominated by large veloc- 
ity dispersion halo members, suggests that the absolute 
number of halo stars in each field is roughly constant. 
However, in field SOUTH, galfast gives an unreasonably 
large halo contribution relative to this expectation, so 
in that case we simply use the Besangon value. For 
field NORTH18, the halo median/ dispersion were fit by 
examining the top quartile of stars in S/N (62 stars, 
S/N> 96), which displayed a relatively narrow peak at 
[Fe/H] = -1. 57. For field NORTH , we used SDSS proper 
motions from iMunn et all ((2004) to isolate 7 halo stars 
with the cut \v ra d\ > 100 km s _1 and /!/ > 220 km s — 1 . 
This yielded a median metallicity of —1.80. For field 
SOUTH, the kinematic data were of lower quality , and 
we adopt [Fc/H] = -1.6 from iCarollo et al.l (|2010f) . For 
fields NORTH, SOUTH the dispersion was calculated by 
adding the errors due to Poisson noise and stellar param- 
eter scatter in quadrature with a 0.3 dex intrinsic halo 
[Fe/H] dispersion (108). 

3.6.3. Fitting Monoceros Parameters 

The normalization of the Monoceros distribution is 
fixed by requiring that the components sum to the to- 
tal number of targets. This leaves two free parameters 
characterizing the Monoceros MDF: its median and dis- 
persion. All components were assumed to be intrinsi- 
cally gaussian. To account for detection inefficiency at 
very low [Fe/H], we modulated the distributions by a 
linear ramp defined to be zero at (and below) the mini- 
mum detected [Fe/H] and unity at (and above) the maxi- 
mum rejected [Fe/H] . For each field, we fit the Monoceros 
[Fe/H] median /j. mon and dispersion a mon by maximizing 
the likelihood: 



L= - V^05io{P( Fe/H }i;/J, m on,<?mon)} (1) 

n 

i=l 

Where n is the total number of reliable [Fe/H] measure- 
ments in each field and P([Fe/H]) is the normalized sum 
of the halo, thick disk and trial Monoceros modulated 
gaussian distributions. Table lists the fitting assump- 
tions and results, which agree with expectations based 
on visual inspection. 

4. CALCIUM ABUNDANCE OF THE MONOCEROS 
STREAM 

We have also measured [Ca/H] using the CallK line 
for program stars sufficiently metal poor to avoid satura- 
tion effects ([Fc/H] < —0.5). Our [Ca/H] determination 
follows a procedure very similar to that described for 
[Fe/H] in i j3.2l For each star, we adopt the same stel- 
lar parameters as in the [Fe/H] fits, but now input the 
measured [Fe/H] as the Kurucz atmosphere metallicity. 
To determine each star's [Ca/H] value, we minimize the 
same x 2 statistic, this time for 3922A < A < 3945A. The 
line list for Ca II K syntheses is provided in Table [3J 

4.1. Continuum Normalization & Calcium Calibration 

For the Ca fitting region, we used a local continuum 
normalization procedure, iteratively fitting a second or- 
der polynomial to the region 3845A < A < 4045A. We 
found that this continuum essentially matched the val- 
ues/slopes of the continuum on the two "sideband" re- 
gions blueward/redward of the Call H/K absorption. 
This fits the continuum shape well, but underestimates 
the normalization by neglecting absorption which causes 
the continuum to be rarely if ever realized. For this rea- 
son, these continuum fits caused [Ca/H] to be system- 
atically underestimated, with the effect becoming more 
pronounced with increasing [Fe/H]. 

We recalibrate the continuum by requiring that our 
thick disk sample exhibit a trend of [Ca/Fe] vs. [Fe/H] 
matching that from the literature. The expected trend 
was computed usin g 150 Milky Way field stars from 
IGratton et al.l (|2003|) . We found good agreement in slope 
could be achieved by mulitplying the initial continua by 
a factor of (1 + 6), with 6 increasing linearly from zero at 
low metallicity, [Fe/H] M13 = -1.55, to 13% at [Fe/H] = 
—0.50. Agreement in offset was achieved with subsequent 
addition of 0.03 dex to all [Ca/Fe]. Applying this proce- 
dure to our 189 M13 spectra yield ed [Ca/Fe] = 0.25 
in agr eement with the estimate of ICohen fc Melendea 
(f2005h . who find [Ca/Fe] = 0.19 ± 0.07 from high- 
resolution spectroscopy of 25 cluster members. 

We apply the local continuum fit and subsequent recal- 
ibration to Monoceros candidates with — 1.1 < [Fe/H] < 
—0.5, where we had sufficient disk stars to verify the 
recalibration procedure. Figure 0] overplots these Mono- 
ceros results with the literature trend and recalibrated 
disk, where the lines are constructed by taking the me- 
dian of [Ca/Fe] in bins 0.2 dex wide in [Fe/H]. Also 
plotted are linear fits to recent [Ti/Fe] abundances 
of 70 Sagittarius and 21 Monoceros stream members 
(jChou et al.ll2010aHH ). The a-enhancements [Ti/Fe] and 
[Ca/Fe] are expected to display similar trends, as both 
reflect the relative heavy element contributions of Type 
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Fig. 3. — (top row, bottom left) Decompositions of the [Fe/H] apec distributions for each of our three Monoceros target fields. In each 
case, the green line represents the assumed halo contribution, the purple line represents the assumed thick disk contribution, the red line 
represents the best fit Monoceros contribution, and the cyan line gives the sum of these three components. The black histogram shows 
the observed [Fe/H] spe c distribution in bins 0.3 dex wide, while the blue line is a curve which counts the number of objects within a 0.3 
dex wide [Fe/H] bin, as the bin center is shifted incrementally by 0.01 dex at a time, (bottom right) [Fe/H] spe c distributions of thick disk 
targets in Galactic north (thick black line) vs. south (thin black line) displaying the observed 0.22 dex metallicity asymmetry. 



II (a-elements+Fe) vs. Type la supernovae (Fe only), 
and hence probe the star-formation history. dSph galax- 
ies typically show deficiencies in a-elements relative to 
the Milky Way, owing to a slower star formation rate in 
which fewer Type II explosions occur before the onset of 
Type la explosions at >1 Gyr. 

5. DISCUSSION 

5.1. Monoceros Abundances 

As the Monoceros members comprise ~2/3 of targets 
in each field, it can be seen by visual inspection that 
Monoceros has a median [Fe/H] ~ — 1, distinct from both 
the thick disk and halo. Indeed, our maximum likelihood 
fits yield [Fe/H] ~ —1.0 for each subsample, showing no 
trend with I (see Tabled]). 108 analyzed ~ 11, 000 MSTO 
stars in a similar spatial region to that of our targets 
(13 < i?cc/kpc < 16, 3 < Z/kpc < 4, 170° < I < 190°, 
see their Figure 18); our inferred Monoceros medians 
-1.09 < [Fe/H] < -0.98 are consistent with the 108 
photometric estimate of [Fe/H] = —0.95 to within the 
systematic uncertainty of ~0.13 dex on our spectroscopic 
abundances. It should be noted that we have used the 
same g — r temperature calibration as 108. A more strin- 



gent test of SDSS photometric metallicities would bepro- 
vided by a study deriving T e g independently. See §5.3l for 
a detailed comparison of photometric vs. spectroscopic 
metallicities. 

108 also inferred a small intrinsic Monoceros [Fe/H] 
dispersion of only 0.15 dex. For each field, we simulated 
the distributions of [Fe/H] arising from a monometal- 
lic population at [Fc/H] = \x mon , given our S/N dis- 
tribution, T e ff distribution, continuum placement errors 
and expected stellar parameter scatter. We found RMS 
measurement-induced scatter of 0.21 dex, 0.23 dex, and 
0.16 dex for fields NORTH, SOUTH, NORTH18, imply- 
ing intrinsic scatters of 0.21 dex, 0.22 dex, 0.10 dex re- 
spectively. Applying the same methodology to M13 stars 
of comparable S/N (>29), we find a measurement error 
of 0.36 dex, and observed scatter of 0.42 dex, yielding 
an intrinsic M13 scatter of 0.20 dex. Globular clusters 
gene rally have RMS iron abundance dispersion < 0.05 
dex (jCarretta et al.l [20091 ). However, with only a sin- 
gle exposure, the M13 sample is subject to additional 
sources of error (e.g., cosmic rays) which we have not 
modeled, perhaps reconciling this difference. If we are 
indeed slightly underestimating the observation-induced 
scatter, a very conservative interpretation would regard 
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TABLE 2 

[Fe/H] Distribution Fitting Assumptions/Results 



[Fe/H] 

Field N S i ars N haio N T d N mo „ JiTD cttd (dexj Jih^To ^halo (dex) /^ m0 n %on (dex) L maa; 

NORTH 183 40.8 23.2 119.2 -0.89 0.35 -1.80 0.42 -1.09 0.30 -0.21 

SOUTH 131 33.3 13.0 84.7 -0.67 0.35 -1.60 0.42 -0.98 0.33 -0.23 

NORTH18 247 28.7 56.4 161.8 -0.89 0.35 -1.57 0.25 -1.03 0.19 -0.10 



the intrinsic Monoceros dispersions we obtained as upper 
limits. 

For fields NORTH, NORTH18, a single gaussian com- 
ponent is sufficient to describe the Monoceros MDF. 
For field SOUTH, there appears to be an excess at 
[Fe/H] ~ — f .2. This could arise from a multimodality in 
the Monoceros MDF or, more likely, indicates contami- 
nation by the metal weak thick disk (MWTD, see §6.2j) . 
Inspection of the CMD in field SOUTH also suggests 
multiple main sequence turnoffs aside from that of the 
disk. Given these uncertainties in population modeling, 
the inferred gradient in Monoceros median metallicity 
with b (+0.11 dex from b = +25° to b = -25°) appears 
not to be significant. 

In two of three fields, the Monoceros [Ca/Fe] dips be- 
low the Galactic trend as [Fe/H] increases from —1.1 to 
-0.5 (see Figure 3). In fields NORTH and SOUTH, this 
[Ca/Fe] trend is qualitatively similar to the recently mea- 
sured [Ti/ Fe] trends of Monoc eros and Sagittarius stream 
members (jCnou et al.l 120 10bl |al). In field NORTH18, the 
observed Monoceros trend appears to exactly match that 
of the Milky Way However, the field NORTH18 decom- 
position in Figure [3] suggests the measured [Fe/H] dis- 
tribution cuts off more sharply than expected at high 
metallicity. If indeed [Fe/H] has been systematically un- 
derestimated for the highest metallicity stars in this field, 
[Ca/Fe] would be artificially inflated. It should be noted 
that our Monoceros sample plotted is a mixture of Mono- 
ceros and disk stars, driving down the apparent contrast 
between the Monoceros and Milky Way trends. 

5.2. Thick Disk Abundances 

In addition to the Monoceros candidates, in fields 
NORTH and SOUTH we measured [Fe/H] for 30 and 
60 thick disk stars respectively. These fields exhibit a 
strong stellar number density asymmetry, with the south- 
ern field showing an excess of more than 50% compared 
to its northern counterpart. The asymmetry is easily 
detected in the SDSS data, is also detectable in Pan- 
STARRS1 37r survey data, and cannot be accounted for 
by photometric errors or reddening uncertainties. J uric 
et al. (in preparation) speculate it may be a signature of 
yet another, more nearby, stream, or otherwise a global 
disturbance (a "warp" or a "bend" ) in the thick disk, of 
uncertain origin. 

Figure [3] shows the measured distributions; the median 
and dispersion of [Fe/H] in field SOUTH are -0.65, 0.29 
dex respectively, while for field NORTH, the correspond- 
ing values are —0.87, 0.22 dex. The statistical scatter 
in our [Fc/H] mcasurments ~0.05 dex per star has neg- 
ligible effect on these values. The median \Z\ for the 
NORTH, SOUTH thick disk targets is respectively 1.60 
kpc and 1.79 kpc, using the photometric parallax rela- 
tion of 108 equations (A6) and ( A7). According t o the 
gradient in disk metallicity from iBond et al.l (|2010D this 



corresponds to an expectation, given perfect disk sym- 
metry, that SOUTH would have [Fe/H] 0.01 dex lower 
than NORTH. The asymmetry is not simply reflecting 
different scale heights. Although the adopted systematic 
error on our metallicities is nominally 0.13 dex, this rep- 
resents an offset with respect to the true abundances; our 
internal consistency should be much better, especially 
for these two thick disk samples which have very similar 
colors and S/N, as well as relatively strong Fe features. 
Reconciling the median metallicities of NORTH/SOUTH 
thick disk sample would require that stellar parameters 
such as logg and v t differ greatly between the two pop- 
ulations, which itself would require invoking disk asym- 
metry about the midplane. A Kolmogorov-Smirnov test 
between the two measured [Fe/H] distributions yields a 
consistency probablility of <2.5 x 10~ 3 , even excluding 
the added 0.01 dex contrast from the \Z\ gradient in 
[Fe/H]. 

The dispersion in [Fe/H] for field SOUTH is 0.22 dex, 
versus 0.29 dex for field NORTH and 0.35 dex for the 
thick disk (averaged over many directions) according to 
iBensbv et "all (|2007t ). A disrupted dSph origin for the 
overdensity would naturally explain the lower [Fe/H] dis- 
persion in field SOUTH, but also incorrectly predict an 
[Fe/H] lower than that of the northern disk for the south- 
ern stellar overdensity. 

Thus, we have spectroscopically identified an asymme- 
try in the thick disk towards the Galactic anticenter, at 
-Rgc ~ 12 kpc and \Z\ w 1.7 kpc, which we expect to 
be confirmed and characterized further via stellar counts 
analysis. 

5.3. Comparison with Photometric Metallicities 

With ugriz photometry avail able, we can comp ute 
[Fe/R] phot from equation (Al) of IBond et all (|2010f ) for 
the entire sample. jFe/HL j,, t is reddening-sensitive. 
ISchlaflv fe Finkbeinerl (|2010h have recalibrated SFD, 
finding the need for a substantial, ^45%, global change 
in E(u — g), but only a 2% change in E(g — r). Thus, 
our T e ff are negligibly affected at a <10K level. But for 
fixed g - r, and near [Fe/B] phot = -1, [Fc/B] phot has 
a str ong u — q gradient, >0 - 04 de x per 0.01 mag u — g. 
The ISchlaflv fe Finkbeinerl (|2010l ) recalibration implies 
A(u — g) >0.06 in portions of our sa mple, a major effect . 
However, the [Fe/Hjp^ot estimator of IBond et al.l (|2010f ) 
is tied to SFD reddenings, and we therefore calculate 
[Fe/H]pft t using (u — g)sFD] it may be necessary to re- 
vise the SDSS [Fe/H] p fto( relation using updated E{u — g) 
values if these photometric metallicities are to be inde- 
pendent of dust column. 

Furthermore, single-pointing samples of [Fe/H] p / lot and 
[Fe/H] spec are sensitive to local reddening in differing 
ways ([Fe/H] S p ec is a function of g — r reddening through 
T c ff ). Near [Fe/H] = —1, using conventional redden- 
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Fig, 4. — (l e ftmos t three panels) [Ca/Fe] as a function of [Fc/H] for Monoceros fields. The green line is the Galactic trend from 
IGratton et al. (2003), the blue line represents our thick disk sample, and the red line represents Monoceros candidates. All lines are 
constructed by taking the median of individual measurements in bins 0.2 dex wide in [Fe/H] centered about the abscissa. Open circles 
are measurements for individual stars, and the crosses indicate typi cal single-star errorba rs at the low metallicity end, [Fe/H] = —1.1. 
(rightmost panel) The green line shows the [Ti/Fe] trend inferred from Gratton ct al. (2003) in a fashion identical to that used to construct 
the green [Ca/Fe] lines at left. The dot-dashed a nd dashed red line s give linear fits to [Ti/Fe] measurements of 21 Monoceros members 
and 70 Sagittarius stream members, respectively (Chou ct al. 2010a EJ). Two of three Monoceros fields (NORTH, SOUTH) show [Ca/Fe] 
underabundance trends relative to the Milky Way similar to those displayed by [Ti/Fe] in Sgr and Mon. 



ing laws, overestimating A r leads to an underestimate 
of [Fe/H] p;iot , but to an overestimate of [Fe/H] spec for 
fixed absorption EW. Any local (degree scale) error in 
A r drives our photometric and spectroscopic metallici- 
ties apart. 

Nevertheless, we find a highly significant correlation 
between [Fe/H] p ;, ot and [Fe/H] spec . The correlation is 
strongest for the highest S/N spectra, corresponding to 
the brightest objects with least noisy ugr photometry 
and best measured [Fe/H] spec . M13, closely approximat- 
ing a monometallic population, shows only a marginal 
correlation, in accordance with the expectation that the 
observed spread in [Fe/H] values arises predominantly 
from random scatter in the photometric/spectroscopic 
measurements. 

Figure [5] overplots the distributions of photomet- 
ric/spectroscopic metallicities for each of our subsam- 
ples. In general, the median photometric and spectro- 
scopic metallicities agree at the level of ^0.15 dex, which 
can be reconciled by the estimated systematic uncertain- 
ties of 0.13 dex on [Fe/H] spec and 0.1 dex on [Fe/H] p /j i 
(108). Except for the M13 sample, the [Fe/H] spec distri- 
butions are narrower than their [Fe/Hjp/jot counterparts, 
by a median of 0.07 dex in dispersion. Of course, much 
of the total dispersion is intrinsic. The superiority of 
M13 [Fe/Hjphot owes to the targets being relatively bright 
(r ~ 18), yielding good photometry, yet receiving only a 
single, short Hectospec exposure. 

5.4. Radial Velocities 

As a byproduct of our abundance extraction procedure, 
we have heliocentric radial velocities in hand for our en- 
tire sample. Typical statistical errors are 10-20 km s _1 , 
depending on S/N. We gauged our systematic error using 
M13 as a calibration sample, since t his globular c luster 
has a known v rad = -244.2 km s _1 (|Harrisl[l996} ). We 
find for our 189 M13 spectra a median of v ra d = —241.1 
km s _1 , and therefore have applied a correction of —3.1 
km s _1 to all radial velocity measurements. 



5.4.1. Monoceros Candidates 

lYannv et al.1 ()2003[ ) found that Monoceros is kinemati- 
cally cold, consistent with the expectation for a disrupted 
dSph galaxy, measuring radial velocity dispersions o~ r = 
22-30 km s _1 . For stars with ^ mo „ — a mon < [Fe/H] 
< Mmcm + o-mon in fields NORTH, SOUTH, NORTH18, 
we find cr r =35.7, 31.3, 26.3 km s^ 1 after correcting for 
median measurement uncertainties of 17, 24, 11 km s _1 
respectively. For comparison, the galfast (Besancon) pre- 
dictions for a r of thick disk stars in these fields is 57.7 
(38.4), 56.6 (36.3), 54.3 (35.8) km s -1 . All dispersions 
were calculated as half of the difference between the 
84th and 16th percentile v ra d values. In all cases, the 
metallicity-selected Monoceros sample is kinematically 
cold relative to the thick disk model, and appreciably 
so if we to are adopt the galfast predictio n. Our range 
of ay is slightly higher than that found by lYannv et al.l 
(2003), but this may be due to contamination of the sam- 
ple by disk and halo stars, both of which would tend to 
increase ay. 

We can also compare the median v rac i for the 
metallici ty-selected Monoceros s amples to the predic- 
tions of Pena rrubia et al.1 (|2005l ). Figure [6] overplots 
measured radial velocities with t he th ree prograde 
orbit models of iPefiarrubia et al] (|2005f ). Figure [6] 
also plots I vs. b for o ur Monoceros targets and 
the IPefiarrubia et all (|2005f ) orbits, showing that field 
SOUTH is on a different "wrap" of the stream than are 
NORTH, NORTH18. There is indeed qualitative agree- 
ment between all of the models and our measured me- 
dian v rad for NORTH, SOUTH, NORTH18 of 7 ± 2.9, 
49 ± 4.5, —20 ± 2.4 km s _1 , trending upwards with in- 
creasing I. Quantitative agreement is only achieved at 
the ~10 km s _1 level, with no single model particularly 
favored, as shown by the histograms of the metallicity- 
selected Monoceros v ra d measurements in each field (Fig- 
ure E}. The median Monoceros v rac i values generally dif- 
fer substantially from the galfast (Besangon) thick disk 
v ra d predictions of 30 (22), 54 (59), 5 (6) km s _1 for fields 
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Fig. 5. — Photometric (magenta) vs. spectroscopic (blue) [Fe/H] distributions. The distributions generally agree at the level of ~0.15 
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NORTH, SOUTH, NORTH18 respectively. In summary, 
the measured radial velocities of our metallicity-selected 
Monoceros subsamples are kinematically cold relative to 
the thick disk, have different median values than ex- 
pected for thick disk stars, and genera l ly agr ee with the 
model predictions of iPenarrubia et ail ([2005D . 

5.4.2. Thick Disk Targets 

Since, as discussed in N5.21 we find that the thick disk 
SOUTH sample is chemically distinct from its northern 
counterpart, it is of interest to compare the measured 
thick disk radial velocities with ga^/asi/Besancon pre- 
dictions. For both the NORTH and SOUTH thick disk 
samples, we find rather poor agreement in terms of me- 
dian v ra d, with measured values of —8, 31 km s _1 ver- 
sus galfast (Besancon) predicted values of 21 (18), 52 
(49) km s _1 . The measured dispersions a r for NORTH, 
SOUTH are 42, 25 km s _1 , versus galfast (Besangon) 
predicted values of 46 (47), 46 (49) km s _1 . The mea- 
sured velocity dispersion in field NORTH agrees with the 
thick disk predictions, whereas the SOUTH targets are 
much more kinematically cold than expected. Thus, the 
thick disk kinematics also display a north/south asym- 
metry, with a low velocity dispersion for the southern 
overdensity suggestive of a coherent stream. 

6. CONCLUSIONS 

6.1. What is the Monoceros Iron Abundance? 

We have presented the first ever spectroscopic Mono- 
ceros MDF based directly on Fe absorption line mea- 
surements, finding [Fe/H] = —1.0 ± 0.1. Our re- 
sults confirm the photometric met allicity analyses of 108 
([F e/H] = -0-9 5) and [§Isaret_aL| (|2^TlT ). who also used 
the lBond et al.l (|2010l ) calibration to derive results consis- 
tent with [Fe/H] 1.0 towards (/,6) = (232°,26°). How- 
ever, these earlier results are vulnerable to systematic 



problems (e.g., u band extinction uncertainties) which 
are largely circumvented by our spectroscopic analysis. 

While C03 and Y03 inferred [Fe/H] for Monoceros 
spectroscopically, neither measured Fe absorption di- 
rectly. Instead, both relied upon indices calibrated to 
Mg and Ca features. Our [Ca/Fe] results suggest that 
such calibrations of [Fe/H] to standard a-element trends 
may not be justified. Further, in the case of Y03, whose 
technique measured Ca II K EW, our analysis shows that 
Ca and Fe absorption in this region are highly degener- 
ate, implying large uncertainties for the Y03 procedure. 
There does not appear to be a simple means of bringing 
about agreement between our result [Fe/H] = —1 and 
that of Y03, given the similar location of their targets at 
(Z,&)=(198° -27°), d Q = 13 kpc and our field SOUTH 
at (Z,b)=(203°,-24°), d Q = 12 kpc. To whatever extent 
the present results disagree with those of C03 and Y03, 
our results should take precedence, as we have measured 
[Fe/H] directly, analyzing spectra with depth-of-exposure 
~7x that of Y03 and S/N comparable to that of C03. 
Our sample size also eclipses that of C03 by a factor 
>10x and that of Y03 by a factor ^1.5x. 

Assuming the Monoceros progenitor contained stellar 
populations of varying age, we can still reconcile our 
result with that of the C03 M-giant sample, [Fe/H] = 
—0.4 ± 0.3. M-giants preferentially trace yo unger, higher 
metal licity populations than MSTO stars. iCarlin et al.l 
(|2011| ) have noted that this bias reconciles their MSTO 
Sgr stream [Fe/H] = —1.15 with literature M-giant val- 
ues of [Fe/H] ~ —0.6. The Monoceros case appears di- 
rectly analogous given our conclusion that Monoceros has 
MSTO [Fe/H] = -1.0. This reasoning can also bring 
about agreement between our [Fe/H] = — 1 and the me- 
dian [Fe/H] = —0 .71 of the 21 kinem atically-selected M- 
giants studied by iChou et~aLl (|2010bt ). 

Since we have measured [Fe/H] directly from the Fe 
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Fig. 6. — Adaptation of Pcnarrubia et al. (2005) Figure 7. (bottom) Prograde orbit models (prol, pro2, pro3) of Monoceros debris 
position in the anticenter direction. The gap at I ~ 195° signifies that field SOUTH corresponds to a different orbital "wrap" than 
NORTH, NORTH18. (top) Model-predicted radial velocity as a function of galactic longitude, with each point representing a star in our 
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radial velocities by field. Capped, solid vertical lines represent 
measured v ra d medians and their standard errors. Corresponding 
(uncapped) vertical lines represent the model predictions at each 
field location, with the same linetype legend as in Figure [6] 



absorption lines of MSTO stars (a population with rel- 
atively little metallicity bias) and with S/N and sam- 
ple size comparable to or better than previous spec- 
troscopic efforts, we recommend adoption of our value 
[Fe/H] = —1.0 ± 0.1 as the metallicity of the Monoceros 
stream. 



6.2. The Metal Weak Thick Disk 

iCarollo et all (|2010f ) have speculated that the MWTD 
and Monoceros may share a common origin. This no- 
tion is base d in part on t h e coin cidence between the 
preliminary IWilhelm et al.l (|2005f) estimate of Mono- 
ceros' metallicity, [Fe/H] = —1.37, and their inferred 
([Fc/H]) mwtd = —1.3. Our measurements argue 
against this MWTD-Monoceros connection, as we find a 
~0.3 dex offset in [Fe/H] between these populations. All 
of our systematics are characteristically ^0.05-0.10 dex, 
and crudely speaking we rule out a median [Fe/H] = — 1.3 
at the 2-3a sys t level. Decreasing our median metallic- 
ity to [Fe/H] spec = — 1.3 would require, for example, an 
unreasonably large temperature scale miscalibration of 
-330K. 

6.3. The Monoceros Progenitor 

Much controversy and speculation has surrounded the 
origin of the Monoceros structure since its discovery. To 
characterize the stream's nature, Y03 used star counts 
to estimate the total stellar mass of Monoceros at ^2 x 
IO'Mq-5 x 10 8 M Q , a range consistent with the content 
of a re latively luminous dwarf galaxy. iPenarrubia et al] 
( 2005) specifically simulated the case of a dSph progen- 
itor, fitting observations well with an initital progeni- 
tor of total mass 3 -9xlO s M P) . Still, detractors (e.g., 
IMomanv et al.l [2006D argue that there is no need for an 
extragalactic origin of Monoceros. How can our obser- 
vations further constrain the nature of the Monoceros 
progenitor? 

The trends of ([Fe/H]) and its dispersion, cpe/H], 
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Fig. 8. — (left) ([Fe/H]} vs. My, where each circle represents a single dwarf galaxy, with properties quoted from INorris ct al. (2010j) and 
referefences therein. Frrorbars are only shown when the uncertainty on {[Fe/H]) exceeds 0.1 dex. The dashed line shows our linear fit to the 
eight plotted galaxies brighter than My=— 8. The blue rectangle represents the allowed region of parameter space, combining this computed 
trend with our measured MDF, and con strains the Monoceros progenitor to have —13.2 < My < —12.3. Alternatively, we can use the 
Monoceros progenitor mass estimated by Pcnarrubia ct al. (2005) and assume the Sgr mass-to-light ratio to plot our Monoceros MDF (red 
hatched rectangle). Clearly, Monoceros follows the general dwarf galaxy trend of i ncreasing ([Fe/H]) with increasing luminosity, (right) 
CT [Fe/H] vs - My. Circles are individual dwarf galaxies with properties drawn from lNorris et all l|2010l ) and references therein. Errorbars 
are only shown when the uncertainty on o"[ Fe / H ] exceeds 0.05 dex. Again, we have plotted Monoceros as the red hatched rectangle by 
assuming the Pcnarrubia ct al. (2005) mass and Sgr mass-to-light ratio, and as a blue hatched rectangle using the {[Fe/H]} constraint on 
My. Monoceros also follows the dwarf galaxy trend of decreasing <T[ Fe/ / H j with increasing luminosity. 



with dwarf galaxy luminosity have been studied by 
INorris et alj (|2010), and are reproduced in Figure [51 
Dwarf galaxy ([Fe/H]} increases monotonically with lu- 
minosity; by fitting a linear trend in ([Fe/H]) vs. My 
for —2.1 < ([Fc/H]} < —0.9, we can translate our 
MDF measurement ([Fe/H]) w — 1.0 ± 0.1 into a con- 
straint on the Monoceros progenitor, — 13.2<ATy< — 12.3, 
assuming the progenitor was a dwarf galaxy (see the 
dashed line and blue rectangle in Figure [5]). Using 
a typical stel lar mass-to-light ratio fo r dwarf galaxies 
(M/L v )*«3 (IChilingarian et all [20ll . ([Fe/H]) = -1 
implies 3.3 x 1O 7 M of stellar mass in the Monoceros 
progenitor. Provided the fraction of the Monoceros pro- 
genitor wd2ich_lia£_be«ristnpped is similar to that of Sgr 
(>2/3. K fc Majewskfl20l0t ). we conclude that our es- 
timate of the stellar mass in the Monoceros stream based 
on chemistry and the assumption of a dwarf galaxy ori- 
gin is consistent with the independent Y03 estimate from 
observed star counts. Having constrained the progenitor 
luminosity, we can also place Monoceros on the Figure 
[5J plot of My vs. <7[ Fe / H ], finding very good qualititave 
agreement between Monoceros (blue rectangle) and other 
bright, low metallicity dispersion dwarf galaxies. 

Instead of deriving My from ([Fc/H]), we can 
check whether a d Sph progenitor with b est-fitting 
mas s determined by IPenarrubia et "ail (|2005t ) conforms 
the INorris et all (|2010T ) MDF trends. Taking the 
Monoceros stream to have the same total mass-to- 
light ra tio as the Sgr core (My =—13.3, M fof =2.5 x 
10 8 M^. iMaiewski et all l200l lLaw fc Maiewskl l2010h . 
the IPenarrubi a et all (|2005l ) mass range translates to 
— 14.7<My< — 13.5. We have combin ed this constraint 
with o ur measured MDF to plot the IPenarrubia et al] 
(2005) Monoceros progenitor (red rectangles) in Figure 
[51 In both cases, My vs. ([Fe/H]} and My vs. crr Fe / H ], 
Monoceros shows qualitative agreement with the dwarf 
galaxy trend. 

Thus, we conclude that our abundance study supports 
a dwarf galaxy origin of Monoceros, in that: (1) The 



MDF median and dispersion both conform to the ex- 
pectations for an appropriately luminous dwarf galaxy 
progenitor, and (2) We have detected [Ca/Fe] deficien- 
cies typical of dwarf galaxy stellar populations in two of 
three fields observed. 

6.4. Future Extensions 

In the future, similar observations/analysis could be 
applied to Monoceros fields that more optimally search 
for I and b abundance gradients in the stream. For exam- 
ple, fixing I and observing two \Z\ values symmetrically 
above/below the Galactic midplane would better isolate 
Monoceros abundance trends with b. Detecting or ruling 
out a gradient toward the midplane in this way would 
significantly constrain models whi ch claim the stream t o 
be a disk "warp" or "flare" (e.g.. IMomanv et alll2006l) . 
Observing an extended range of I values at fixed b would 
better constrain a metallicity gradient along the stream. 
For example, Sgr exhibits such a gradient, but only at the 
level ~2.4xl0~ 3 d ex per degree along the debris trail 
(|Keller et al.l [20101. To probe a gradient at this level 
with spectroscopy comparable to that reported here, ob- 
servations would need to span ^100° in I, a factor of four 
increase relative to our current analysis, 178° < I < 203°. 
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TABLE 3 
Fe A5269, CaIIK Line Lists 



A (A) 


Fe A5269 




A (A) 


CaIIK 


l°g9/ 


A (A) 


CaIIK cont'd 


logff/ 


Species 


E.P. (eV) 


Species 


E.P. (cV) 


Species 


E.P. (cV) 


5268.6 


Till 


2.59 


-1.96 


3922.1 


Fel 


3.29 


-2.44 


3933.3 


Scl 


2.3 


-1.29 


5268.9 


Ti I 


3.32 


-1.74 


3922.7 


Fe I 


2.99 


-2.0 


3933.4 


Scl 


0.02 


-0.65 


5269.5 


Fel 


0.86 


-1.32 


3922.9 


Fel 


0.05 


-1.65 


3933.6 


Fe I 


3.07 


-1.16 


5269.9 


Til 


1.87 


-1.74 


3923.0 


Scl 


1.99 


-2.69 


3933.6 


Fel 


3.27 


-2.05 


5270.3 


Cal 


2.52 


0.16 


3923.0 


Fel 


3.25 


-2.27 


3933.7 


Call 


0.0 


0.11 


5270.4 


Fel 


1.6 


-1.34 


3923.5 


Sell 


0.32 


-2.41 


3933.9 


Til 


2.29 


-2.32 


5271.6 


Til 


2.77 


-0.87 


3924.5 


Til 


0.02 


-0.94 


3934.2 


Til 


0.05 


-2.14 


5272.0 


CrI 


3.44 


-0.42 


3925.2 


Fel 


3.29 


-1.4 


3935.3 


Fe I 


2.85 


-1.87 


5273.2 


Fel 


3.29 


-0.99 


3925.6 


Fel 


2.83 


-1.03 


3935.8 


Fe I 


2.83 


-0.88 


5273.4 


Fel 


2.48 


-2.16 


3925.8 


Scl 


1.99 


-2.27 


3935.9 


Fe I 


3.27 


-2.13 


5273.4 


CrI 


3.45 


-0.7 


3925.9 


Fel 


2.86 


-0.94 


3936.3 


Til 


2.48 


-2.73 


5274.4 


Til 


3.09 


-3.03 


3926.0 


Fel 


3.24 


-0.93 


3936.5 


Scl 


2.32 


-1.6 


5274.6 


Til 


2.42 


-2.97 


3926.3 


Til 


2.58 


-1.91 


3936.8 


Fe I 


3.25 


-1.94 


5275.2 


CrI 


3.37 


-0.35 


3927.7 


Scl 


1.99 


-1.36 


3937.0 


Til 


2.29 


-2.38 


5275.3 


CrI 


2.88 


-0.28 


3927.8 


Scl 


1.99 


-1.69 


3937.3 


Fe I 


2.69 


-1.46 


5275.8 


CrI 


2.88 


-0.05 


3927.9 


Fel 


0.11 


-1.52 


3937.6 


Til 


2.3 


-2.34 


5276.0 


Fell 


3.19 


-2.21 


3927.9 


Fel 


2.83 


-2.27 


3937.7 


Fe I 


3.02 


-1.82 


5276.1 


CrI 


2.88 


-0.1 


3928.1 


Fel 


3.21 


-0.93 


3937.9 


Til 


2.17 


-2.28 


5277.0 


Til 


3.33 


-1.67 


3928.8 


Til 


2.31 


-1.71 


3938.0 


Til 


2.27 


-2.08 


5278.2 


Scl 


3.05 


-1.74 


3929.0 


Til 


2.04 


-1.5 


3938.6 


Scl 


2.3 


-1.6 


5278.2 


Til 


2.34 


-2.5 


3929.0 


Til 


2.3 


-1.85 


3938.6 


Scl 


2.34 


-1.72 


5278.3 


Scl 


3.05 


-3.05 


3929.1 


Fel 


2.76 


-1.88 


3939.7 


Til 


2.31 


-2.47 


5279.7 


Fel 


3.3 


-3.44 


3929.2 


Fel 


3.25 


-1.34 


3940.0 


Scl 


2.61 


-2.12 


5280.3 


CrI 


3.36 


-0.73 


3929.9 


Til 


0.0 


-1.06 


3940.1 


Til 


2.5 


-2.96 


5281.8 


Fel 


3.03 


-0.83 


3930.2 


Scl 


1.99 


-1.4 


3940.9 


Fel 


0.96 


-2.6 


5282.4 


Til 


1.05 


-1.3 


3930.2 


Til 


1.5 


-1.23 


3941.3 


Fel 


3.27 


-1.01 


5283.4 


Til 


1.87 


-0.49 


3930.3 


Fel 


0.09 


-1.49 


3941.4 


Scl 


2.0 


-1.5 


5283.6 


Fel 


3.24 


-0.43 


3930.5 


Fel 


3.21 


-2.52 


3942.4 


Fel 


2.99 


-2.07 


5284.1 


Fell 


2.89 


-3.2 


3930.9 


Fel 


2.45 


-2.86 


3942.4 


Fel 


2.85 


-0.95 


5284.4 


Til 


1.04 


-2.43 


3931.1 


Til 


2.29 


-2.04 


3942.8 


Fel 


3.27 


-1.31 


5285.0 


Scl 


2.5 


-1.42 


3931.1 


Fel 


3.27 


-1.14 


3943.3 


Fel 


2.2 


-2.35 


5285.6 


CrI 


3.36 


-1.13 


3931.3 


Fel 


3.24 


-1.9 


3943.7 


Scl 


2.0 


-2.74 


5285.8 


Scl 


2.5 


0.38 


3931.6 


Til 


2.32 


-1.63 


3944.3 


Scl 


2.61 


-2.36 










3932.0 


Scl 


1.99 


-1.34 
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-1.74 
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